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Themophotovoltaic (TPV) conversion of IR radiation emanating from a radioisotope heat source is under 
consideration for deep space exploration. Ideally, for radiator temperatures of interest, the TPV cell must convert 
efficiently photons in the 0.4-0. 7 eV spectral range. Best experimental data for single junction cells are obtained 
for lattice-mismatched 0.55 eV InGaAs based devices. It was suggested, that a tandem InGaAs based TPV cell 
made by monolithically combining two or more lattice mismatched InGaAs subcells on InP would result in a 
sizeable efficiency improvement. However, from a practical standpoint the implementation of more than two 
subcells with lattice mismatch systems will require extremely thick graded layers (defect filtering systems) to 
accommodate the lattice mismatch between the sub-cells and could detrimentally affect the recycling of the 
unused IR energy to the emitter. A buffer structure, consisting of various InPAs layers, is incorporated to 
accommodate the lattice mismatch between the high and low bandgap subcells. There are evidences that the 
presence of the buffer structure may generate defects, which could extend down to the underlying InGaAs layer. 

The unusual large band gap lowering observed in GaAsi_ x N x with low nitrogen fraction [1] has sparked a new 
interest in the development of dilute nitrogen containing lll-V semiconductors for long-wavelength optoelectronic 
devices (e.g. IR lasers, detector, solar cells) [2-7]. Lattice matched Gai_ y lnyNxAs-|. x on InP has recently been 
investigated for the potential use in the mid-infrared device applications [8], and it could be a strong candidate for 
the applications in TPV devices. This novel quaternary alloy allows the tuning of the band gap from 1.42 eV to 
below 1 eV on GaAs and band gap as low as 0.6eV when strained to InP, but it has its own limitations. To 
achieve such a low band gap using the quaternary Gai_ y ln y N x Asi_ x , either it needs to be strained on InP, which 
creates further complications due to the creation of defects and short life of the device or to introduce high content 
of indium, which again is found problematic due to the difficulties in diluting nitrogen in the presence of high 
indium [9]. An availability of material of proper band gap and lattice matching on InP are important issues for the 
development of TPV devices to perform better. To address those issues, recently we have shown that by 
adjusting the thickness of individual sublayers and the nitrogen composition, strain balanced GaAsi_ x N x /lnAs-i_ y N y 
superlattice can be designed to be both lattice matched to InP and have an effective bandgap in the desirable 0.4- 
0.7eV range [10,11]. Theoretically the already reduced band gap of GaAs-|. x N x , due to the nitrogen effects, can be 
further reduced by subjecting it to a biaxial tensile strain, for example, by fabricating pseudomorphically strained 
layers on commonly available InP substrates. While such an approach in principle could allow access to smaller 
band gap (longer wavelength), only a few atomic monolayers of the material can be grown due to the large lattice 
mismatch between GaAs-i_ x N x and InP (~3.8-4.8 % for x<0.05, 300K). This limitation can be avoided using the 
principle of strain balancing [12], by introducing the alternating layers of lnAs-i_ y N y with opposite strain (-2.4-3. 1% 
for x<0.05, 300K) in combination with GaAs-i.xNx. Therefore, even an infinite pseudomorphically strained 
superlattice thickness can be realized from a sequence of GaAs-i_ x N x and lnAs-i_ y N y layers if the thickness of each 
layer is kept below the threshold for its lattice relaxation. Figure-1 shows the band edge variation (red curve) in 20 
period of GaAs^NVInAs^Nx superlattice lattice matched to InP. As we have discussed earlier, the quaternary 
alloy GalnAsN could be a viable material for the application in IR devices, but as seen under the band 
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anticrossing model approximation [13] the band gap of InGaAsN lattice matched on InP (blue curve) is much 
higher than that of the superlattice for the same nitrogen content. Hence the given superlattice band gap can be 
achieved for much lower nitrogen concentration in comparison to that in InGaAsN quaternary. This is of great 
technological advantage, due to the fact that the higher nitrogen concentration is hard to incorporate in the alloys 
and it also deteriorates the quality of the alloy. 



Figure 1. Evolution of the band edge energy gaps at 1 0K in a 20 period GaAs 1 _ x N x (3ML)/lnAsi_ x N x SL as a 
function of the nitrogen composition in the SL (red curve, after Ref.10). Note for each different nitrogen data point 
the thickness of the InAs^Nx sub-layers (4-9MLs) is adjusted to satisfy the lattice matching of the SL to InP (001). 
For comparison purpose band edge of quaternary GalnAsN is also shown (blue curve). 


Figure 2 b) shows the schematic of the insertion of the superlattice in the / region of the conventional p-l-n diode 
shown in Figure 2a), with a schematic of the superlattice in figure 2c). Also shown in figure 2c) is the formation of 
miniband structures in conduction and valance bands, which determines the effective band gap of the structure. 
The effect of insertion of such a superlattice-like alloy within the intrinsic region of a 0.74 eV InGaAs p-l-n diode 
was previously evaluated and it was shown that such a single junction device exhibits a photovoltaic response 
comparable to its lattice-mismatched 0.55 eV-InGaAs counterpart as shown in figure 3 [after Ref. 14]. In this work 
we have extended the approach to multi-junctions devices. Here three or more subcells with different effective 
bandgaps for the superlattice region are monolithically series connected. Maximum power output and 
performance of double, triple and quadruple junction TPV cells are evaluated as a function of the superlattice/cell 
design at a given black body emitter temperature. The study stresses the potential of the proposed approach for a 
significant enhancement of TPV converter performance. 

Following the treatment of the Olson et.al. [15], we can make following simplifying assumption for the well 
behaved p-l-n diode: 1) transparent zero resistance tunnel-junction interconnects, 2) no reflection losses, 3) no 
series resistance losses, 4) junctions collects every photon absorbed and 5) l-V curves are described by the ideal 
(n=1) equation. The short circuit current density (j sc ) of i th subcell is determined by the quantum efficiency of the 
subcell, Q(A) and by the spectrum of the light incident on the that cell, (/{A) as, 
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Figure2. Schematic of the p-l-n diode, defining different parameters, a). Conventional p-l-n device, bj. p-l-n 
device, with /-region filled with a GaAsN/lnAsN superlattice, c) Schematic of the superlattice and the minibands. 
Vertical arrow in c) determines the distance between the miniband edges to give the effective band gap of the 
superlattice structure. 

J‘c =q\ Q Qi(X)fa)tU ( 1 ) 

Based on the assumptions we have made, quantum efficiency can be assumed simply dependant on the total 
thickness, x of the device, as 

Q, W = 1 - exp[- «, {X)x i ] (2) 

because the fraction exp[-a(/l)x] of the incident light is transmitted through the cell instead of being absorbed. 
For the photons with wavelength greater than the band gap, a(A)~ 0, and hence exp[-a(/l)x] = 1 . The incident 

radiation flux (f{X) on the top cell is simply the flux i/> BB (A) of the black body radiation hitting the top surface of the 
cell. On the other hand, the photon flux hitting the bottom or lower cells is filtered by the top cells, so the bottom 
cell only sees an incident spectrum reduced by the factor of exp[-a 1 (/l)x 1 ] of t/> BB (A), where ai(A) and Xi are the 
absorption coefficient and the thickness of the top or the first cell. In general the short circuit current of i th cell in 
the multijunciton cell comprised of m subcells can be given as, 


J'sc = (l - exp[- «, [X)x l ])exp 


in 


<t> BB (A)dA 


i=2 


( 3 ) 
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For example, in three-cell system, the above relation for top (/= 1 ) and bottom (/=3) cell becomes, 



Energy, eV 


Figure 3. (After Ref. 14). Increased photocurrent response due to the addition of superlattice in / region of p-l-n 
device. For comparison purpose effects of InGaAs lattice matched and lattice mismatched to InP are also shown. 

Ac = < : /£ l (i - ex p[ - «i Uk \)t>cH AW- (4) 

and 

Ac = c l\' t] (l - exp[- a 3 (/l)x 3 ])exp[- (a, (A)x, + a 2 (X)x 2 )}j) BB {l )dA (5) 

respectively. We can assume that last cell is infinitely thick to absorb all the photons, giving, for example, J S c of 
last cell in three-junction cell as 

Ac = ex Ph («i Uk + «2 U ) x 2 (Add (6) 

Classical ideal photodiode J-V equation can be written as, 
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where the dark current density J 0 is given as, 

Jo=Jo, base + Jo, intrinsic+J 0. emitter 

with different terms given as, 
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Here, L is the diffusion length, S is the surface recombination speed, n, is the intrinsic carrier concentration in the 
intrinsic region, D is the diffusion coefficient and the distance parameters x and w are defined in Fig. 2. Radiative 
recombination coefficient B rad is carrier density independent property of the material and is given as [1 6], 


B..„ = 


' rad 


n 2 h 3 c 2 


IXJ 

j a(E)z 


E/kT 


E 2 dE 


( 12 ) 


where, n s is refractive index of the material, h is the Planck’s constant, c is the speed of light. From equation-7, 
the voltage across the i th subcell is given as, 
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and the voltage across the tandem device comprised of m subcells can be given as, 


m U T m 

V to Jj)=£v i (j)=— 
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Total power produced by this device can be given as, 


( 14 ) 
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Tandem Current (A) 


Figure 4. Transcendental expression of equation 17 plotted as function of current. Transcendental function of 
equation=0 determines the optimal current through the system. 

Definition of the efficiency of the TPV system is generally very vague and depends on the device, radiation 
source, structure, converter and source separation and many other considerations [17,18]. One specific example 
of efficiency is radiative heat conversion efficiency as described by Mahorter et.al. [19]. In this study instead of 
efficiency we have focused on the power output density, an approach adopted by many authors , for example in 
ref. [20], to determine the device performance, since it depicts the performance of the converter itself, which is the 
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main intention of this work. We consider the blackbody radiator operating at a temperature of 1350K and the 
tandem temperature as 300K. In the first approach the band gap of first cell is fixed to 0.74eV, which is the band 
gap of the InGaAs lattice matched to InP to find the optimized current through the system giving the maximum 
power output by varying the band gap of the second subcell. Current through the first cell is 0.64A/cm 2 and by 
solving the transcendental equation 17, as shown in figure 4, optimal current for the two subcell tandem is found 
to be 0.618 A/cm 2 when the summation T(J) goes to zero. Once the optimal current for the given set of tandem is 
found, equation-15 can be used to determine the output power of the tandem. For a fixed first subcell bandgap, 
second subcell bandgap can be varied to find the maximum power output of the tandem, giving the corresponding 
second subcell bandgap as the optimal band gap of the double junction tandem. Following the procedure, we can 
find the band gaps of the second (0.625eV), third (0.535eV) and fourth (0.46eV) subcells to achieve maximum 
power output in double, triple and quadruple junction cells respectively. Figure 5 shows the result of the procedure 
discussed above for the double, triple and quadruple junction devices, giving the output power density variation 
with the change of the band gap of the last subcell. Solid curves show the optimal powers for the subcell 
consisting bulk like lattice mismatched InGaAs in the conventional p-l-n structure, on the other hand, the broken 
curves are the calculations using the superlattice in the / region of the p-l-n structure. We can notice the better 
output performance using the superlattice, in addition to the lattice matching on the InP substrate. Blue and red 
dots respectively shows the results of the theoretical and experimental work performed by 
Wanlass et al [20] and Wilt et al [21]. 



Last cell band gap (eV) 


Figure 5. Power output for double (blue), triple (green) and quadruple (red) junction cells. Maximum power on 
each curve gives the optimal band gap for second (0.61 5eV), third (0.525eV) and fourth (0.46eV) subcells for 
fixed first (0.74eV) subcell. Solid curves are for the subcells with lattice mismatched InGaAs materials and the 
broken curves are calculated including the superlattice materials in the / region of the p-l-n device. 

If the first band gap is not fixed, optimal band gaps of the first two subcells also can be found by searching all the 
possible combination of the band gaps of the first subcell and the second subcell, giving the iso-power surfaces or 
iso-power curves as shown in figures 6 and 7 for the lattice mismatched InGaAs devices and in figure 8 for the 
superlattce device. It can be seen that the optimal band gaps for first and second subcells comes out to be at 
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0.64eV/0.46eV for the lattice mismatched InGaAs devices and 0.62eV/0.41eV for the superlattice devices 
respectively. It is important to note that the power output density of the tandem device made of superlattice (fig 8) 
is much higher than that of the lattice mismatched InGaAs (fig 7). 



0.35 0.4 0.45 0.5 0.55 

Band gap of second cell (eV) 


Figure 6, Iso-power surface for the band gap variation of top two subcells in double junction tandem consisting of 
lattice mismatched InGaAs. Colored scale shows the power in W/cm 2 . It can be seen that maximum output power 
can be reached with slight variation for different band gap combinations. 

As we have discussed that the lower band gaps needed for the subcells of multijunction device can be obtained 
from an lattice mismatched InGaAs layers on InP, with consequence of poor device performance due to the strain 
related problems. Replacing the / region with unstrained lattice matched GaAsN/InGaAs SL would avoid the 
presence of strain in the device on the one hand but on the other hand would reduce the active absorbing layer 
thickness, which is limited to only the / region of the subcell, since p and n region would be an unstrained InGaAs 
with fixed band gap of 0.74eV. One consequence of thinner active absorbing region would be lesser current 
density and lesser power output. Since the absorption coefficient of the GaAsN/lnAsN SL is higher than that of 
bulk InGaAs [14], / region filled with a superlattice will perform better than the / region of InGaAs, but since the 
overall active absorption region would be decreased, there will be a trade off between the strain effects on the 
device and the maximum power output. Since for the many applications, longevity of the device is an issue, it is 
important to have stable devices, even compromising the overall device output yield. With a use of relatively 
thicker superlattice we can gain the power density output equivalent or even better in comparison to the lattice 
mismatched devices. For a quadruple cell, the power density output reached more than 0.8W/cm 2 for the 
superlattice thickness of 0.2pm for the second and third subcell and ~1 pm for the fourth subcell, which further can 
be increased by optimizing the thickness and other physical properties of the individual cell. Main reulsts of the 
paper can be summarized in the figure 9. This figure shows the variation of the output power density of the 
double, triple and quadruple tandem device with and without the use of a superlattice. The optimization of both 
the band gaps for double junction tandem with and without the superlatices are shown as green and black dots 
respectively. Hence higher output power can be reached using the optimized double band gaps with the use of 
the superlattice. 
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Figure 7. Iso-power curves representing combination of band gaps of first and second subcells made of lattice 
mismatched InGaAs to produce a given output power. Maximum power output can be reached for band gaps of 
0.64 and 0.46eV for first and second subcells respectively. 
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Figure 8. Iso-power curves representing combination of band gaps of first and second subcells made of 
superlatice to produce a given output power. Maximum power output can be reached for band gaps of -0.62 and 
0.41 eV for first and second subcells respectively. 
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Number of cells in Tandem 


Figure 9. Output power density vs. the number of tandem cells for the devices with and without the use of lattice 
matched superlattice. Green and black dots respectively show the output power for the double junction tandem 
with and without the superlattice, in which both the band gaps were optimized for the highest power output. 

In conclusion we have proposed the use of strain balanced superlattice structure in TPV device to achieve variant 
bandgaps between 0.4-0. 7eV as needed for the type of radaitive source used. We have shown that better or 
comparable performance of device can be achieved by replacing the / region of the subcell by strain balanced 
GaAsN/lnAsN superlattice lattice matched to InP. Due to the lattice matching structures, it is technologically 
feasible to implement three or four junction cell on InP based devices, which otherwise in the case of lattice 
mismatching devices is hard to implement due to the requirement of thick graded layer between the cells. We also 
deduced the optimized band gaps for two-cell tandem device to produce optimal power. 
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